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times of pulses. This effect can also arise in a nonlinear material. We propose a model in which 
fluctuations of a background electric field, such as that produced by a squeezed photon state, can 
cause fluctuations in the effective lightcone for probe pulses. This leads to a variation in flight 
times analogous to that in quantum gravity. We make some numerical estimates which suggest 
that the effect might be large enough to be observable. 
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I. INTRODUCTION 



A quantum theory of gravity is expected to predict quantum fluctuations of the space- 
time geometry, and hence of the hghtcone. This effect was recognized long ago by early 
workers on quantum gravity [1-3]. In particular, Pauli [1] once suggested that lightcone 
fluctuations might remove the ultraviolet divergences of quantum fleld theory by smearing 
the singular behavior of Green's functions. This hope has not yet been realized, but still 
remains a possibility, given that a complete quantum theory of gravity has not yet been 
found. Several proposals for the small scale structure of spacetime near the Planck scale 
have been made, including that of "spacetime foam" [4, 5]. For a recent discussion of work 
on this topic, including observational bounds, see Ref. [6]. The possibility of an energy- 
dependent refractive index in one model was discussed in Ref. [7]. Even in the absence of 
a full theory, it is possible to discuss lightcone fluctuations in weak field quantum gravity. 
When linearized perturbations of a fixed background are quantized, it becomes possible to 
treat smearing of Green's functions and the associated lightcone fluctuations [8-10]. These 
are the active fluctuations, coming from the dynamical degrees of freedom of gravity itself. 
There are also passive fluctuations, which are driven by the fluctuations of the stress tensor 
of matter flelds [11]. 

One expects quantum gravity effects to be very small, except in extreme conditions such 
as the early universe or near small black holes. For this reason, it is of interest to seek analog 
models in condensed matter systems, where similar effects might be exhibited but be much 
larger. The use of analog models for black hole evaporation has been a very active area of 
research, and is reviewed extensively in Ref. [12]. The effects of spatial randomness in a 
material have recently been treated as an analog model for lightcone fluctuations [13, 14]. 
In the present paper, we wish to propose a different model based on nonlinear optics, where 
fluctuations of the effective dielectric function lead to a form of hghtcone fluctuations. 

There has long been interest in the parallels between light propagation in gravitational 
fields and in dielectric media. Gordon [15] showed that a gravitational field may be used 
to mimic the effects of a dielectric. Conversely, it is possible to define a dielectric medium 
which mimics the effect of a gravitational field [16]. In the latter case, the medium must 
have the property that the magnetic permeability is equal in magnitude to the dielectric 
permittivity in Gaussian units. This is not realized by any known material. Our interest 
is not in attempting to reproduce the detailed effects of gravity in a medium, but rather in 
the fact that a realistic material alters the speed of light in the material. When this speed 
can fluctuate, a form of lightcone fluctuations arises. 

The outhne of the paper is as follows: In Sect. II, we review selected aspects of lightcone 
fluctuations in linearized quantum gravity. The crucial results from nonlinear optics which 
are needed will be discussed in Sect. III. In Sect. IV we introduce our analog model, and 
derive a key result for the variance in photon flight times in the presence of the lightcone 
fluctuations. Some numerical estimates are given in Sect. V, where we discuss the possibility 
of a laboratory experiment to look for the lightcone fluctuation effect. Our results are 
summarized and discussed in Sect. VI. We will employ SI units, except in Sect. II, where 
h = c = 1 units are used. 
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II. LIGHTCONE FLUCTUATIONS IN LINEARIZED QUANTUM GRAVITY 



In this section, we review selected aspects of lightcone fluctuations in quantum gravity. 
First, let us recall the situation in classical general relativity theory, where the lightcone 
is a dynamical object determined by the spacetime geometry. A small perturbation of the 
geometry alters the lightcone. If we start with a background of Minkowski spacetime, and 
add small perturbations, the apparent speed of light as measured by the flat background 
metric can either increase or decrease. In the case of a weak field Schwarzschild metric, the 
change is a decrease, leading to the well-known Shapiro time delay [17]. The measurement 
of the time delay of radar signals passing near the sun is one of the better experimental tests 
of general relativity [18]. Similarly, a Schwarzschild metric with a negative mass would lead 
to a time advance, or an apparent increase in the speed of light compared to flat spacetime. 
Of course, the local speed of light is constant; it is the apparent speed measured over a finite 
distance which changes. 

Now suppose that the linearized perturbations arc quantized and are subject to quantum 
fluctuations. This will lead to fluctuations in the apparent speed of light and hence of the 
hghtcone. Let the metric be written as 



(1) 



where r]^^ is the Minkowski metric, and hf^^, is the quantized linear perturbation, which now 
becomes the graviton fleld operator. Let a{x,x') be the invariant interval function, deflned 
as one-half of the squared geodesic distance between spacetime points x and x'. In flat 
spacetime, a — ao, where 

ao^kx-x'r. (2) 



Expand a in powers of hn^ as 



a = ao + ai + a2 + 



(3) 



so that (Ti is first order, a2 is second order, etc. In flat spacetime, the retarded Green's 
function for a massless field is a delta function on the forward lightcone: 

e{t - 1' 



Go{x,x') 



(4) 



In a fixed curved spacetime, the form of the retarded Green's function near the lightcone 
has the same functional form, but with ctq replaced by a. If we now suppose that /i^^, and 
hence ai, undergo Gaussian fluctuations around a mean value of zero, then the averaged 
Green's function becomes a Gaussian centered on the mean lightcone [8]: 



{G{x,x'))^ 




exp 



2{af) 



(5) 



In the limit that (af) — )■ 0, the Gaussian in Eq. (5) approaches a delta function, and the 
averaged Green's function approaches the form in Eq. (4). In this case, the mean lightcone 
is that of Minkowski spacetime. If we include the second order term, (72, then Eq. (5) is 
modified to become 



{G{x,x')) 




exp 



(^0 + (^2))' 



(6) 
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Wc notice that these two last equations appear misprinted in Ref. [8]. The effect of the (T2 
contribution is simply to shift the mean lightcone from that of flat spacetime to a curved 
spacetime lightcone. One way to create lighcone fluctuations is with a bath of gravitons in a 
non-classical state, such as a squeezed vacuum state [8]. In this case, (cti) = 0, so to lowest 
order, the lightcone fluctuates around that of flat spacetime. The shift of the mean lightcone 
described by ((T2) reflects the curving of the background spacetime by the gravitons. 

Here (af) and ((J2) are understood to be finite expectation values. If we are interested 
in the effect of gravitons in a non-classical state, then we may take the difference between 
the formal expectation values in the state in question, and in the Minkowski vacuum state. 
More generally, finite expectation values can arise from averaging two-point functions over 
finite regions of spacetime. 

The replacement of a delta function by a Gaussian is the mathematical expression of 
lightcone fluctuations. The physical consequence of the smearing of the hghtcone is that the 
speed of propagation of pulses between a source and a detector becomes statistical. Some 
pulses will travel slower than the mean speed, as measured in the background spacetime, 
but others travel faster, with a Gaussian distribution of flight times centered on the classical 
flight time in flat spacetime. This leads to a mean variation in flight times of 

At^^, (7) 

r 

where r is the distance (measured in the flat background) between the source and the 
detector. This variation may be expressed in terms of the graviton two-point function 

{h^^{x)Ki3{x')) as 

(At)2 = ^ / dtdt' n^rfn'^n^ {h^,{x)Kp{x')) . (8) 

Here is a unit spacelike vector defining the direction between the source and detector, 
and the integrations are along the mean lightcone. This expression may be shown to be 
gauge invariant provided that the metric perturbation h^^y is localized between the source 
and the detector [19]. Equation (8) may also be derived from the Riemann tensor correlation 
function [20]. Integrals of the Riemann tensor correlation function may also be used to obtain 
other physical effects of spacetime geometry fiuctuations, such as luminosity variations [21] 
and spectral line broadening and angular blurring [22]. 

It is important to note that At is the ensemble averaged variation in flight times, not 
necessarily the expected variation in flight times of two successive pulses [9]. If the two 
pulses are sufficiently close in time, they are correlated and tend to have flight times which 
differ by less than At. This can be understood as correlated pulses probing approximately 
the same classical spacetime geometry. A given source of metric fluctuations, such as a bath 
of gravitons, will deflne a characteristic fluctuation time scale of r. Pulses emitted at time 
intervals less than about r are correlated, whereas those emitted with longer separations 
become uncorrelated and have flight time differences of order At. If the metric fluctuations 
are due to a bath of gravitons with characteristic wavelength A^, then pulses are uncorrelated 
when their times of emission differ by more than about A^. 

For the case of a squeezed vacuum state of gravitons, pulses traveling in the same direction 
as the gravitons have a flight time variation of the order of 

AtPi lp\g\fUr , (9) 
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where ^p is the Planck length, and U is the characteristic energy density of the gravitons. 
Clearly, this is a very small effect in the present day universe. The presence of the Planck 
length in Eq. (9) can be traced to the fact that the graviton two-point function is proportional 
to £p, and is hence a universal feature of lightcone fluctuations in quantum gravity. A major 
motivation for analog models is to find systems with larger effects. 



III. WAVE PROPAGATION IN NONLINEAR OPTICS 



A. Wave Equations 

Unlike general relativity, Maxwell's theory is fundamentally a linear theory. However, 
interactions of electromagnetic fields with charges can produce an effective nonlinear wave 
equation. In nonlinear optics, this arises because atoms behave as nonlinear oscillators in 
strong electric fields. For a review, see for example, Ref . [23] . In the absence of free charges 
or currents, electromagnetic waves in a nonlinear material medium is described by the usual 
Maxwell equations: 

V-B = 0, VxE = -— , V-D = 0, VxH=— , 10 

dt dt ^ ^ 

together with the constitutive relations for a nonmagnetic, but polarizable material, 

B = /ioH, D = eoE + P. (11) 

Here E and B stand for the fundamental electric and magnetic fields, respectively, while D 
and H are the corresponding induced fields. The relation between the polarization vector 
P and E is nonlinear, and can be expanded in a power series as 

= eo {x^^E, + x^IEjE, + x%E,E,Ei + •••), (12) 

where sums on repeated indices are understood. The coefficients of the various powers of 
the electric field are the components of the susceptibility tensors. Here Ei stands for the 
total electric field, which could be in part due to the polarization itself and in part due to an 
applied external field. We have assumed that the polarization at a time t depends only on 
the instantaneous value of the electric field. This is a reasonable approximation in a regime 
where dispersion and dissipation can be neglected. 

The wave equation for the electric field follows from Eqs. (10) and (11) and is 

V(V.E)-V=E + -L^E=--1,^P, (13) 

where c = 1/ ^/(/ioeo) is the speed of light in vacuum. We will be interested in cases where 
V • E = 0, so the wave equation becomes 

V2_1^)e=^— P. (14) 
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We will restrict our attention to solutions of this equation which are polarized in the z- 
direction, but propagating in the x-direction, so that 



Ei ^5i^E ^ E{t, x) . 
In this case, V • E = is satisfied. Now we may write 

P.-eo{x^'^E + x^'^E' + /'^E' + ...) , 

where 



y(l) ^ y(2) ^ (2) 

A A.ZZ 1 A A.ZZZ 1 



y(3) ^ y(3) 
A A,zzzz 



(15) 



(16) 



(17) 



Note that in introducing this notation, we are making no assumptions about isotropy. In 
general the susceptibility tensors have other nonzero components which are not equal to 
the components listed above. Rather, the components listed in Eq. (17) are the only ones 
needed for our discussion. 

Here we will consider nonlinear effects through third order, and ignore any higher order 
effects. In this case, the wave equation for E becomes 



^E--—E--—(y^^^E^ + y^^^E^)^0 



where 



V — 



(18) 



(19) 



is the speed of fight in the 2;-direction when only the linear effects are included. Let -E'o(i, x) 
be a solution of Eq. (18). Consider a second solution of the form E ~ Eq -\- Ei, where 
\Ei\ <^ |£'o|, but El varies more rapidly in time than does Eq, so that 



dEn 



dt 



El 



Er 



dEi 



dt 



(20) 



Assuming x^^-'-Ei << 1 and x'^^^^^o-^'i << 1 and using the approximation that Eq is slowly 
varying, Eq. (18) reduces to the following wave equation for Ei, 



d'^Ei 
dx^ 



- ^ (1 + 2ei + 3e2) 



d'^Ei 



0, 



where 



and 



£2 



X 



(2) 



1 + x(^) 



^(3) ^ 

nr Eq {t, x) 



(21) 



(22) 



(23) 



1 + x(^) 

Note that the subscripts on ei and 62 refer to the power of Eq, not the order of the suscep- 
tibility. 
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B. WKB Solutions 



Equation (21) describes waves which encounter a space and time dependent effective 
dielectric function. In the short wavelength limit which we consider, this translates into 
a propagation speed which varies in space and time. This can be made more precise by 
examining WKB solutions. We consider plane wave solutions of the form 

Ei{t,x) ^ Ne'''^'' F(t) , (24) 

where N is a normalization factor. Equation (21) implies that 

-g^+uj'F^O, (25) 

where 

TT, (26) 

l + 2ei(t) + 3e2(t) ' ^ ' 

and oji = kiv is the angular frequency in the limit that ei = £2 = . A WKB solution of 
Eq. (25) is 

F = F{t,x) = (2uj)-^/^ exp -i [ dt' u(t') 

L Jo 

^ {2uj)-^/^ expi^-icuit + iuji dt' ei{t' , x) + i^^- dt' [e2{t' , x) - el{t' , x)]^ .{27) 
The corresponding solution for the electric field is 

Ei{t,x) — Ni exp i^ikix — iujit + iuji J dt' ei{t' , x) + i^-^ J dt' [e2{t' , x) — el{t' , x)]^ , 

(28) 

where A^i = N{2uj) and we are working to second order in Eq. 
This solution describes a wave propagating with phase velocity 

u{t,x) ^ ^===^==== ^ V {l - e,{t,x) - h2{t,x) - el{t,x)]\ . (29) 
^yl + 2el{t,x) + 3e2{t,x) L ^ J 

In order to discuss the flight times of pulses, we need to form wavepacket solutions. So long as 
dispersion may be ignored, these wavepackets will have group velocities also approximately 
given by Eq. (29). Recall that this space and time dependent velocity is determined in part 
by the background field, Eo{t,x). This is analogous to the situation in general relativity 
theory, where the apparent speed of light depends upon the spacetime geometry. Just as in 
gravity, the effective lightcone for the propagation of weak disturbances depends upon the 
background field. 



IV. FLUCTUATIONS OF THE EFFECTIVE LIGHTCONE 



A. Flight Times of Pulses 

Our analog model for quantum lightcone fluctuations arises when the background field 
Eo{t,x) is no long a fixed, classical field, but is allowed to undergo quantum fluctuations. 
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This leads to quantum fluctuations of the speed of light in the medium and of the effective 
lightcone. The operational meaning of these lightcone fluctuations is a variation in flight 
times of pulses between a source and a detector. If the spatial separation between the source 
and detector is r, then the flight time for nonzero background field is 



r= / — -- — -dx. (30) 

Jq u{t{x),x) 

Here x = v {t — to), where t = to is the time of emission of the pulse. Hence the local velocity 
u{t{x),x) is taken to be evaluated on the unperturbed path of the pulse. When |ei| <^ 1 
and 1 62 1 <^ 1, this may be written as 



T 



-+- r 



ei + ^(3e2-e?) 



dx , (31) 



with ei and 62 understood to be evaluated ai t — to + x/v. 

Now we suppose that Eo{t,x), and hence ei and €2, undergo quantum fluctuations. The 
mean flight time is 

^^^""v^Uo {^'"^ ^ ~ ''0 "^"^ ■ ^^^^ 

and the variance in flight time is 

{Atf = (T^) - (r)2 . (33) 

We will evaluate this quantity to second order in Eq, which is analogous to second order in 
hf,,, in quantum gravity. This means that we retain terms of order or 62, but drop higher 
order terms such as 6162. The result is 

(Aty = — / dxi dx2[{ei(ti,X2)ei(t2,X2)) - {ei{ti,X2)) {ei(t2,X2))] , (34) 
Jo Jo 

where = to + Xi/v. Note that although a term involving (3e2 — ) appears in (T^) and in 
(T)^, it cancels from (At)^. This is analogous to the situation discussed in Sect. II, where 
(a2) modifies the background, and hence the mean flight time, but it is {af) which describes 
the lightcone fluctuations. 

Now we assume that ei fluctuates around a mean value of zero, so that 

(ei)=0, (35) 

then the last term in Eq. (34) vanishes. This expression may also be written as a double 
integral in time: 

pr/v pr/v 

{Atf= dti dt2{ei{h)ei{t2)), (36) 
Jo Jo 

where we drop the explicit reference to the space dependence. Equation (36) is the analog of 
Eq. (8), with the correlation function for ei playing the role of the graviton correlation func- 
tion. Thus a material with nonzero second-order susceptibility, x^'^\ is needed for an analog 
model which reproduces the features of quantum gravity described in Sect. II. Note that if 
we were to work to higher order in Eq, then terms including (^2(^1) ^2(^2)) and (e^(ti) e^(t2)) 
would also appear in (At)^. These terms would model the effects of quantum stress tensor 
fluctuations, or passive quantum gravity effects, as opposed to the active quantum gravity 
fluctuations modeled by Eq. (36). 
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B. Phase Fluctuations 



In addition to flight time fluctuations for wavepackets, there is another related effect 
associated with hghtcone fluctuations. This is the fluctuations in the phase of a plane wave 
solution, such as Eq. (28). Here we will assume Gaussian fluctuations, which is the case 
when dealing with an approximately free quantum fleld. Let 

$ = + ($) (37) 

be the phase of the wave, so that {(p) = 0. Let the fluctuations of (p be described by a 
Gaussian probability distribution of width a centered at = 0, 

P{p) = exp (-^) . (38) 



The ensemble average of any function of (p is deflned by 

/oo 
f{<p)P{<p)d<p. (39) 
-oo 

One may easily verify that 

{e-n = exp (-^ ((^2)) . (40) 

There is an alternative derivation of Eq. (40) in which the ensemble average is an ex- 
pectation value in a "vacuum-like" quantum state so that (/(</?)) = (■0|/(</')|V')- Here 
"vacuum-like" means that there exists a decomposition of (p, treated as a quantum operator, 
into positive and negative frequency parts, (p = p^ + (p~, such that p>^\ip) = 0. Wc may 
write (<^^) = {p}~^ P>~)- This relation may be combined with the Campbell-Baker-Hausdorff 
formula, e"^"*"^ = e^al^'^l e'^e^, for any pair of operators A and B that each commute with 
their commutator [A, B] to obtain Eq. (40). The two derivations are equivalent because the 
"vacuum-like" states are really squeezed vacua, for which the wavefunction in a Schrodinger 
representation is a Gaussian in some variable. 

If we take the average of the electric fleld, as given by Eq. (28), then the result is 

{E^{t, x)) = A{t) e^['=i^-'^i *+^(*)] , (41) 

where 

m = ^ 

^ Jo 
and 

A(t) = (iVi) e-H(^*)'. (43) 

Here (At)^ is given by Eq. (36) with t = r/v. As before, we work to second order in the 
background fleld, Eq, which means second order in ei but flrst order in 62- The effect of 62 
appears only through its expectation value, and takes the form of a phase shift which can be 
interpreted as a shift in frequency of the wave. In second order, the effects of ei also include 
light cone fluctuations, which here show up as a decrease in the amplitude of the averaged 
wave. 



<P{t) = ^ f dt'{e2{f)-el{t')) (42) 
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C. Quantized Electric Field in a Squeezed Vacuum State 

Now we will regard the background field, Eo(t, x), as a quantum field, which may be 
expanded in terms of modes as 

Eo{t, X) = ^[OkABkA hx{t, X) + a^kACkA f kx{t, x)] . (44) 
kA 

Here GkA are polarization vectors for linear polarization A, the a'^kx and ak\ are photon 
creation and annihilation operators, respectively, and the mode functions are plane waves: 

/kA(t,x)=zy^e^('^— *). (45) 

Box normalization in a volume V is assumed, and a; = T;|k|. 

We will take the quantum state of the electromagnetic field to be a multi-mode squeezed 
vacuum state. This is analogous to the bath of gravitons discussed in Sect. 11 and pro- 
duces lightcone fluctuations through a fluctuating effective dielectric function. A single 
mode squeezed vacuum state, \0, is obtained by acting on the vacuum with the squeeze 
operator [24] 

^(C)=exp|i[ra^-C(a^)l}, (46) 

where ( = re^^ is a complex parameter, and and a are the creation and annihilation 
operators for the excited mode. In this state, one has that 

{a^a) = sinh^ r , (47) 

and 

(a^) = -e*'' cosh r sinh r . (48) 

We now wish to construct a multi-mode squeezed vacuum state by allowing the squeeze 
parameter C to depend upon the mode, but assuming that 

(akAOk'A') = (a'^kAOk'A') = (49) 

unless k = k' and A = A'. This insures that the electric field correlation function will be 
a sum of contributions from the individual modes with no cross terms. As in Sect. Ill A, 
we wish to assume that the excited modes are all at least approximately polarized in the 
z-direction, so that 

ekA ~ z , (50) 

which picks out a specific polarization and allows us to omit the mode label A. 

Next we compute an electric field correlation function as an expectation value of a product 
of 2;-components of electric field operators. We are interested in the effects of excited modes, 
and take the product to be normal ordered and symmetrized to write 

{Eoit,K) Eo{t',^')) = i(: E^{t,^) E^,{t',^') : + : E^,{t',^') E^,{t,K) :) (51) 

^ _L Re V cu I (ak^Ok) e^[k-(x-x')-(*-01 _ (^^2^ g.[k.(x+x')-.(t+t')] | . 
eol/ ^ I J 
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When each of the excited modes is a squeezed vacuum state, this becomes 

h 



{EQ{t, x) EQ(t', x')) — — — ^^cj sinhrk {sinhrk cos[k • (x — x') — uj{t — t')] 

k 



+ coshrk cos[k ■ (x + x') - a;(t + t') + r/k]} • (52) 
The mean squared electric field is the coincidence limit of this expression 

{E^it, x)) = — — V] w sinhrk {sinhrk + coshrk cos[2k • x — 2ujt + ?^k]} • (53) 

Note that this quantity can be negative, which is an example of a subvacuum effect, similar 
to negative energy density in quantum field theory [25]. 

The phase shift, 0, defined in Eq. (42), can be expressed as a time integral of the squared 
electric field as 



0(t) 



3 uji 



2(l+x"l) 



2 ^ rt 



I dt' {Elif ,^)) . (54) 

^0 



1 + X«. 

Similarly, the mean squared time delay, given by Eq. (36), may be written as 

{^tf = (j^) [ jl "^'^^ ■ 

Here wc understand x'(t') and x"(i") to be on the mean path of the light ray. Wc assume that 
the excited modes are all approximately moving in the +a;-direction, so that x(t) ^ vtH. 
However, in the limit of large quantization volume V, there are necessarily still a large 
number of excited modes, and we may write 

fE-(^/'^''=-(^*'A.A«. (56) 



k 



Here Ak <^ k is the bandwidth in wavcnumber of the excited modes, and ^ 1 is their 
angular spread around the +a;-direction. As in Sect. HI A, we assume that all of the relevant 
modes have a phase velocity of v = cu/k, so that 

k • X fa . (57) 

This causes the arguments of the trigonometric functions in Eqs. (52) and (53) to be ap- 
proximately constant. Here the electric field correlation function and squared electric field 
become time- independent, and given by 

h 

(EqU, x) EqU', x')) = (Eq) = — — uj k'^ Ak A9 sinhr(sinhr + cosh r cos?]) . (58) 

Here we have used Eq. (56) and assumed that the squeeze parameter is approximately the 
same for all excited modes, so the mode label on r and rj has been dropped. 

Note that the energy density U in empty space associated with the excited modes is 

U^eo {El) . (59) 
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It will be convenient to express our final results in terms of U . The phase shift becomes 



U 3t 



X 



(3) _ 



i + xW. 



(60) 



Because oc t, it can be interpreted as a shift in frequency of the original probe field of 



Sou 



3U 



26o(l + x(^): 



X 



(3) _ 



1 + x(^) 



(61) 



More interesting for our purposes is the variance in fiight times of pulses, which becomes 



{Atf 



X 



(2) 



U 



(62) 



This result is the analog of Eq. (9), and tells us that the root-mean-squared fiuctuation in 
fiight time can grow linearly with increasing path length. Thus the lightcone fiuctuations in 
a nonlinear material with nonzero x^'^^ can potentially become large. The correlation time 
in this case is expected to be of the order of the period of the excited modes in the squeezed 
state, which form the background field. Thus pulses separated by more than this period 
should have fiight time variations of order At. 



V. SOME NUMERICAL ESTIMATES 



Here we wish to give an estimate of the magnitude of the lightcone fiuctuation effect in 
nonlinear materials and an assessment of the possibility of observing it in an experiment. 
First, let the ratio of the root-mean-squared fiuctuation in fiight time, A^rms — \/ (At)^, to 
the total fiight time be denoted by 5, so that 



5 = %^ = W--^— (63) 




This may be expressed as 



_ 3.36 X 10-6 I U 

l + X^ 10-11 m/V y lJ/m3 • 

We need to make estimates of each of the parameters which appear in Eq. (64). The first 
order susceptibility, x''^^ is typically of order or less than unity at optical frequencies, so the 
1/(1 -l-X^^'') factor will be of order one. The second order susceptibility, x^"^^ is niore variable, 
and can only be nonzero for materials which do not possess spatial inversion symmetry (i.e., 
are noncentrosymmetric) . Such materials typically exhibit a significant degree of anisotropy, 
SO the different components of the tensor xl/k '^^^Y considerably in magnitude. We need a 
material for which a diagonal component {i = j = k) of this tensor is nonzero, which we 
can take to be xizz = X*"^"*- An example of such a material is cadmium sulfide (CdS), for 
which x^^^ ~ 4 X 10-11 m/V. (See Fig. 22 in Ref. [26], noting [23] that xi^l = 2^33.) Other 
materials, such as gallium arsenide, possess somewhat higher values of the second order 
susceptibility, but have only off-diagonal components of Xi% being nonzero [23] . 
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Next we need an estimate of the realizable energy density in a squeezed vacuum state. 
Experiments producing squeezed light have achieved 10 db of squeezing [27], meaning 10 = 
lOlogiQ(U/Uy), where Uy is the vacuum energy density, so that U ^ lOUy for each squeezed 
mode, or about five photons per mode. The number of modes excited in a finite quantization 
volume V is the product of V with the right hand side of Eq. (56). Thus 

where \ = 27i/k is the mean wavelength of the excited modes, which will typically be in the 
infrared or optical ranges. This can be written as 

J /1000nm\^ fAk 



To be consistent with our earher analysis, we need both Ak/k and to be small. If, 
for example, we have Ak/k ~ ^ 10~^, then fractional flight time variations of order 
6 ~ 10^^ would seem to be possible using infrared light. Although this is small, it might be 
observable with accurate flight time measurements. 

Recall from Eq. (43), that the growth of (At)^ with flight time causes the amplitude of 
a plane wave solution to decay exponentially in time. This is a feature of dephasing, where 
phase fluctuations lead, via Eq. (40), to a decaying amplitude. This is also in principle 
an observable characteristic of lightcone fluctuations. However in an experiment, it may 
be very difficult to distinguish this effect from the exponential decay due to other forms of 
absorption of the wave. 

Although gravity is dispersionless, real dielectric materials exhibit frequency dependence 
in their susceptibilities, which we have ignored. In addition, we have assumed that the 
frequency of the probe field is higher than that of the background field. In general, the two 
fields will experience different values of the first and second order susceptibilities, which we 
have not considered. However, many materials with nonzero x^'^\ including CdS, exhibit a 
X^^^(c<;) which is relatively independent of frequency in the infrared and the red end of the 
visible spectrum [26, 28]. The second order susceptibility starts to vary more strongly with 
frequency at higher frequencies. Thus our analysis should give an approximate description of 
lightcone fiuctuation effects in the red and infrared parts of the spectrum. There has been 
recent progress on designing quantum dot hetereostructures with specified second order 
susceptibilities [29]. This raises the possibility of engineering materials with the properties 
needed for this type of experiment. 



VI. SUMMARY AND DISCUSSION 

In this paper, we have argued that nonlinear optical effects provide an analog model for 
quantum lightcone fluctuations in quantum gravity. The basic idea is that a fluctuating 
background electric fleld in a nonlinear material produces a fluctuating effective dielectric 
function for a probe fleld, leading to fluctuations in the effective speed of light. This phe- 
nomenon is not unique to nonlinear optics, but will arise in any system with a nonlinear 
wave equation. Sound wave in a fluid, for example, satisfy the nonlinear Westerwelt equa- 
tion. Thus a fluctuating background sound wave can induce fluctuation of the "soundcone" , 
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and hence leading to a similar analog model. However, nonlinear optics seems to be a more 
promising arena for precision measurements. 

Materials with a nonzero second order susceptibility, x^"^^ > provide the best analog models 
for the effects in linearized quantum gravity. The variance of the flight time of a pulse, 
which can be taken to be the signature of lightcone fluctuations, is proportional to an 
integral of an electric field correlation function in the analog model, while the same quantity 
is an integral of the graviton correlation function in quantum gravity. We considered a 
background electric field created by photons in a multimode squeezed vacuum state, which 
induces quantum fiuctuations in the effective dielectric function experienced by probe pulses. 
Our numerical estimates suggest that fractional fiight time variations of order 10~^ might 
be achievable and potentially measurable. An alternative to photons in squeezed vacuum 
states might be a stochastic bath of electromagnetic waves, for which {Eq) = 0. This might 
allow even larger effects. 

Given that quantum gravity itself is generally remote from experimental tests, further 
study of analog modes such as that proposed in this paper seems to be a promising area for 
further enquiry. 
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